Compressibility Consideration in the Boundary of a Strongly Collapsing Bubble 
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Equations of radial motion of a gas bubble in a compressible viscous liquid have been modified 
to account for compressibility at the bubble boundary. A new bubble boundary equation has been 
derived, including a new term resulted from liquid compressibility. The influence of this term has 
been numerically investigated using isothermal-adiabatic model for the gas inside the bubble. The 
results clearly indicate that, at the end of the collapse, the new term has a very significant role 
and its consideration dramatically changes the bubble characteristics. Moreover, the more intense 
the collapse is, the more significant the effect of the new term is. Also, it has been reasoned out 
that, the influence of the new term will be established even when the effects of mass (water vapor) 
exchange, chemical reactions, and gas dynamics inside the bubble are taken into account. 

PACS numbers: 47.55.Bx, 43.25.Yw, 43.25.+y, 78.60.Mq 



I. INTRODUCTION 

The problem of non-linear radial oscillations of a gas 
bubble in a liquid, when it experiences a high ampli- 
tude spherical sound field, is an old challenging problem. 
Several complications are present in this problem aris- 
ing from the effects of heat conduction, mass diffusion, 
compressibility, chemical reactions, surface tension, and 
viscosity. Many authors have studied different aspects of 
this matter, but a rather complete description has not 
been presented yet. 

The radial dynamics of a bubble in an incompress- 
ible liquid is described by the well-known incompress- 
ible Rayleigh-Plesset equation 0, Q . The extension of 
this equation to the bubble motion in a compressible liq- 
uid has been studied by many authors; for example Her- 
ring Trilling Q , Gilmore Q , Keller and Kolodner [j| , 
Flynn 0, Lastman and Wentzell Lofstedt et al. 
and Nigmatulin et al. ff(|. In these works, several differ- 
ent forms of the compressible bubble dynamics equations 
have been presented, using different approximation for 
consideration of the compressibility effects. On the other 
hand, the influences of heat conduction and mass diffu- 
sion in the bubble motion have been reported by Hickling 
If], Fujikawa and Akumatsu and Yasui |13j . 

In a generalized approach, Prosperetti and Lezzi 0] 
used a singular-perturbation method of the bubble-wall 
Mach number and derived a one-parameter family of 
equations describing the bubble motion in the first order 
approximation of compressibility. This family of equa- 
tions are written as: 



--(Pi-p 




(i) 



where, R, C, Po, P a , and p are bubble radius, liquid 
sound speed, ambient pressure, driving pressure, and 
density of the liquid, respectively. Also, r\ is an arbi- 
trary parameter. Equation must be supplemented by 
a boundary condition equation at the bubble interface to 
relate the liquid pressure, Pi, to the gas pressure inside 
the bubble. Like all previous authors, Prosperetti and 
Lezzi 01 used the following incompressible equation for 
this purpose: 
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where, P g , p, and a are gas pressure at the bubble inter- 
face, liquid first viscosity coefficient, and surface tension, 
respectively. As Prosperetti and Lezzi showed ^4|, dif- 
ferent existing forms of the bubble dynamics equations 
belong to this single parameter family of equations, cor- 
responding to different values of t\. Especially, r\ = I and 
r] = correspond to Keller and Miksis equation and 
Herring- Trilling equation 0, 3 , respectively. 

Two specific approximations have been considered in 
the derivation of Eq'ns. iff} and 10). First, Eq'n. QJ 
has been derived from the Euler equation, in which the 
liquid viscosity effects have been neglected. However, 
these effects are not important in the usual applications 
of the bubble dynamics equations. It can be shown [f6| 
that these effects are significant when p > pCRo, where 
Rq is the bubble ambient radius. Therefore, for micron 
size and larger bubbles the elimination of these effects is 
thoroughly justified. 

The second approximation, which is more important, 
is the incompressibility assumption of the liquid and 
the gas at the bubble interface, which has been used 
in the derivation of Eq'n. J2J). Note that, all of the 
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effects of the liquid compressibility in the work of Pros- 
peretti and Lezzi, as well as in all previous works, have 
been resulted from the liquid motion around the bub- 
ble, but not from the bubble boundary condition equa- 
tion. In fact, all previous authors, on one hand take 
into account the compressibility of the liquid motion 
around the bubble, but on the other hand neglect its 
consideration at the bubble interface. Although Eq'n. 
J2J has widely been used in all old and recent works 
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applicability for the bubble motion needs to be clarified 
(especially at the end of the collapse, where the bubble 
motion is significantly compressible). 

In this paper, we have modified Eq'n. J2J considering 
the effects of compressibility at the bubble interface for 
both the liquid and the gas. The modified equation has 
new terms resulted from the effects of two coefficients of 
viscosity of the liquid and the gas. This work shows that 
the liquid compressibility effects in the bubble dynamics 
originate from two sources. First, they originate from the 
liquid motion around the bubble, which has already been 
presented in the generalized Eq'n. ©. Second, they are 
initiated through the liquid compressibility consideration 
at the bubble interface as presented in this work. We have 
performed numerical analysis to investigate the influence 
of the modification of Eq'n. J2J. 



II. COMPRESSIBLE BUBBLE BOUNDARY 
EQUATION 

To derive the compressible bubble boundary equation, 
we assume that, the motions of the bubble interface and 
the surrounding liquid are always spherically symmetric. 
The continuity equation and the radial component of the 
stress tensor can be written as: 



T rr (liquid) \ R = T rr (gas) \ R + 2 "^- 



(7) 



Applying Eq'n. JSJ for the gas and the liquid parts of 
Eq'n. J7J leads to: 
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where, p g and X g are the first and the second coefficients 
of viscosity of the gas at the bubble interface, respec- 
tively. Also, A; and A g are the divergence of velocity of 
the liquid and the gas, respectively. Substituting the di- 
vergence of velocity for the liquid and the gas from Eq'n 
© into Eq'n. © yields: 
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where p g is the gas density at the bubble interface. Equa- 
tion © represents the bubble boundary equation con- 
taining all effects of the compressibility and the viscosity 
of both the liquid and the gas. Comparison of Eq'ns. 
(J2J) and © indicates the existence of three new terms in 
Eq'n. © due to the liquid and the gas compressibility 
and viscosity effects. Here, we concentrate on the effects 
of the new term arising from the liquid compressibility. 
Therefore, the gas viscosity effects are neglected due to 
their smallness relative to the liquid viscosity as in previ- 

ous works p, 0, n, n n, n, h, ie m m mm, in, 

Under this circumstance, Eq'n. I© becomes: 
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where, p, u, p, and A are density, velocity, pressure, and 
divergence of the velocity, respectively. Also, A is the sec- 
ond coefficient of viscosity. Inserting du/ dr from Eq'n. 
©, into Eq'n. Q yields: 
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The velocity divergence, A, can be written as: 

1 dp 1 dp 

p dt pc 2 dt ' 
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where, the sound speed, c, is defined as c 2 = dp/ dp. The 
boundary continuity requirement at the bubble interface 



To generalize the argument, we express Eq'ns. © and 
(|ll)fl in dimensionless forms. The dimensionless variables 
are: 

r* _ R p* _ R t* - tC p* - Pl 
K 'To' H ~C' To' Fl 'pC 2 "' 



Po 



Pa 



P* = -AS- P* = AA. (II) 

where, Rq is the ambient radius of the bubble. Substitut- 
ing the dimensionless variables into Eq'ns. © and i|l(J|) 
results the dimensionless equations: 



+ 1)R*) R*R* + | (l - l&V + 



)R* ) R* 2 



1 + -rj)R* )(P*-P*-P*) 



(12) 



3 



p; = p;- 4 



p*R* 
R* 



- 2- 



R' 



dP* 

-( A * + V)^f- (13) 



The quantities a*, A*, and p* are the dimensionless 
form of surface tension, first, and second viscosity coef- 
ficients of the liquid and are defined as: a* = a/pRoC 2 , 
X* = X/pR C, and p* = p/pR C, respectively. These di- 
mensionless quantities are basically the inverse of Weber 
Number and Reynolds Number. 

It should be mentioned that, although the effects of 
compressibility in Eq'n. (|12[1 are in the first order ap- 
proximation, but these effects have been introduced com- 
pletely in Eq'n. . To quantify the effects of our mod- 
ification, we have performed our calculations within the 
framework of Eq'n. 1|12[1 . using r\ = 1. Of course, similar 
results are obtained using the other values of 77 in Eq'n. 
(|12|l . because of the similarity of the order of approxima- 
tion in all equations corresponding to different values of 
77. 
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III. NUMERICAL ANALYSIS 

To quantify the effects of the new term in Eq'n. (|13[1 . 
numerical analysis was carried out for the conditions of 
Single Bubble Sonoluminescence (SBSL) [H El. The 
driving pressure in its dimensionless form was Pa{t) = 
— P* sin (ui*t*), where u* is dimensionless angular fre- 
quency and is defined as ujRq/C. To have a well-posed 
problem, the value of the gas pressure at the bubble inter- 
face, P* , must be specified. It can be determined, in the 
most complete approach, from the numerical solution of 
the conservation equations for the bubble interior along 
with Eq'ns. I|12|l and (| 1 31) simultaneously. Also, the bub- 
ble content changes during the bubble evolution because 
of the chemical reactions at the end of the collapse. In 
addition, the mass exchange and the heat transfer be- 
tween the bubble and the surrounding liquid affect the 
bubble interior. 

During the past ten years, several different approaches 
have been presented to describe the real state of the gas 
and its evolution considering the above-mentioned com- 
plexities. These approaches, which at first assumed in- 
viscid and without consideration of chemical reactions, 
heat transfer, and mass exchange 0,H(J, became grad- 
ually more complex by including dissipating effects of 
radiative transfer (2l| , heat transfer, and viscous gas dy- 
namics!^, [2j| . Recent gas dynamic model of Storey and 
Szeri |2J| accounts for all effects of chemical reactions and 
water vapor evaporation and condensation. 

During the bubble motion, strong spatial inhomo- 
geneities inside the bubble are not remarkably revealed 
unless at the end of the collapse. Therefore, the uni- 
formity assumption for the bubble interior seems to be 
useful and provides many features of the bubble motion 
[TtI Il8| . By this assumption, the gas pressure inside the 
bubble is obtained assuming polytropic evolution for the 



FIG. 1: Time variations of the bubble characteristics; (a) 
radius, (b) gas temperature for the compressible boundary 
condition in one period of the applied pressure field. 

gas bubble as: 

where, Van der Waals equation of state has been used. 
Here, a* is the dimensionless Van der Waals hard core 
radius; a* — a/R . The bubble evolution is assumed 
to be isothermal (7 = 1.0) for the radii larger than Rq, 
when the bubble moves relatively slow. While, adiabatic 
assumption (7 = T, where T is the ratio of the specific 
heats of the bubble interior) is applied for the smaller 
radii, when the bubble experiences rapid changes |17| . 

Hilgenfeldt et al. [2f| extended this model considering 
7 to be a function of the Peclet number of the bubble. Of 
course, the Peclet number changes with the variations of 
the bubble radius and the bubble velocity. They showed 
that, many features of SBSL can be explained by their 
model. However, in practice there is very small differ- 
ence (less than 5%) between the Hilgenfeldt et aL's model 
and the isothermal-adiabatic model |2g. Although, these 
two models can not illustrate production of discontinu- 
ities and shock waves at the end of the collapse, but they 
are really useful for our purpose to investigate the impor- 
tance of the new term during the collapse. It should be 
mentioned that, these models have been used in several 
recent works in sonoluminescence studies 12 7j . 

We have used isothermal-adiabatic model in this pa- 
per. Also, we have argued about the extension of the 
results for more complicated models. Under these cir- 
cumstances, the time variations of the bubble proper- 
ties under the framework of Eq'n l|12|) (for r\ = 1) have 
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FIG. 2: Time variations of the bubble radius according to (a) 
compressible and (b) incompressible cases. The time inter- 
val of this figures is 40 ps around the minimum radius time 



FIG. 3: Time variations of the gas temperature when the bub- 
ble reaches its minimum radius according to (a) compressible 
and (b) incompressible cases. 



been numerically calculated for two cases: (a) compress- 
ible boundary condition (Eq'n. I13fl and (b) incompress- 
ible boundary condition (Eq'n. We used Runge- 
Kutta method for our numerical analysis. The con- 
stants and the parameters used in the calculations were 
set for an air bubble in water at room temperature, 
To = 293.0 K, and atmosphere pressure, Pq = 1.0 aim, 
[S3; p = 998.0 kg/m 3 , C = 1483.0 m/s, p = 1.01 x 
10" 3 kg/ms, a* = 1.0/8.745, a = 0.0728 kgs~ 2 , and 
r = 1.4. The second coefficient of viscosity of water was 
set to be A = 2.23 p. It was derived from the value of 
the bulk viscosity of water at room temperature that, 
Pb = A+ (2p/3) = 2.9p [2^. Also, the angular frequency 
of the deriving pressure was u> = 2tt x 26.5 kHz. 



IV. RESULTS 

The results of our calculations have been illustrated in 
Figs. (1-5) for P a = 1.6 atm and Rq = 2.0 pm. Similar- 
values for these parameters were reported in the results of 
recent experimental works of Ketterling and Apfel [30| . 
Simon et al. [31|, and Vazquez et al. |32j. Figure (1) 
shows the variations of the bubble radius and the gas 
temperature in one period of the applied pressure field 
using the compressible boundary condition. The bub- 
ble motion has incompressible characteristics during a 
period, except for an infinitesimal time interval at the 
end of the collapse. Basically, the new term in Eq'n. 
(|13|l originates from the compressibility of the liquid mo- 
tion. Therefore, it is expected that, its effects not to be 



revealed, until the bubble motion becomes significantly 
compressible. The results of our calculations clearly con- 
firm this point. Except for the end of the collapse, the 
differences between the bubble characteristics resulting 
from the compressible and incompressible boundary con- 
ditions are less than 1.0%. This result thoroughly jus- 
tifies the elimination of the new viscous term in Eq'n. 
(|13fl for all times, but not for the end of the collapse. 
Indeed, in Fig. (1), the results for the bubble proper- 
ties of the two boundary conditions completely coincide. 
However, the maximum bubble temperature in Fig. 1(b) 
has a considerable increase relative to that of the incom- 
pressible case. This discrepancy instigated us for further 
concentration of the bubble properties around the mini- 
mum bubble radius. 

In Figs. (2-4), we have presented the evolution of the 
bubble characteristics using compressible boundary con- 
dition (a), and incompressible boundary condition (b), 
around the minimum radius time. The time interval of 
these figures is 40 ps. Since, the bubble experiences its 
maximum compression at the end of the collapse, the ef- 
fects of the new term in Eq'n. (|13fl are more significant 
during this time interval. 

Figure (2) shows the bubble radius evolution for the 
two boundary conditions. This figure illustrates a con- 
siderable difference between the compressible and incom- 
pressible cases. After the bubble reaches its minimum 
radius, a number of small bouncing oscillations appear 
in the graph of the compressible case, which do not oc- 
cur in that of the incompressible case. The period of 
these oscillations is nearly 3 ps. The details of our calcu- 
lations show that, the times of minimum radius for the 
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FIG. 4: Time variations of the dimensionless gas pressure 
(P a /pC 2 ) around the time of minimum radius according to 
(a) compressible and (b) incompressible boundary conditions. 



two cases are the same. While, the minimum radius for 
the new Eq'n. (|T3|> is 1% less than that of the old Eq'n. 

Figure (3) represents the time variations of the gas 
temperature near the minimum radius time. It evi- 
dently represents that, introducing the new viscous term 
in Eq'n. l|13l) strongly affects the gas temperature evo- 
lution at the end of the collapse. After the minimum 
radius time, the incompressible case illustrates smooth 
behavior for the bubble temperature. While, remarkable 
sharp peaks appear in the temperature evolution using 
the compressible boundary condition. Also, as the bubble 
goes away from the minimum radius, the peaks become 
weaker. Note that, the value of the maximum tempera- 
ture increases considerably (nearly 2.4 times) for the new 
boundary equation. 

Comparison of the gas pressure evolution for the two 
cases is shown in Fig. (4). Similar behaviors as in Fig. 
(3) are observed in this figure for the difference between 
the new and the old boundary equations. Note that, the 
gas pressure is much more sensitive to the presence of the 
new viscous term compared to the gas temperature. In 
fact, the presence of the new term increases the maximum 
pressure more than one order of magnitude. 

In Fig. (5), the time variations of the three different 
pressure terms in Eq'n. (|13|l have been illustrated near 
the minimum radius. These terms are due to the ef- 
fects of surface tension and viscosity, i.e. P a = 2a* /R*, 
Pf, = 4n*R*/R*, and P Xfl = (A* + 2[j,*){dP? /dt* ). This 
figure clearly shows that, at the end of collapse, the col- 
lective effects of the viscous terms are by far greater than 
the surface tension term. Moreover, the new viscous 
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FIG. 5: Time variations of three dimensionless pressure terms 
in equation il l -'it namely; (a) surface tension term: P CT = 
2a* j R* , (b) damping viscous term: P M = 4/i* R* / R* , and 
(c) new viscous term: Pa m = (A* + 2fi*)(dP*/dt*), when the 
bubble reaches its minimum radius. 



term, Pa^, is the dominant term in this time interval. 
Note that, the order of the maximum values of the three 
pressure terms are completely different. P\ M is the high- 
est (up to 10 3 ) and P a is the lowest (less than 10~ 2 ). 
These results emphasize that, the elimination of P\^ is 
not reasonable, when the bubble evolves near the mini- 
mum radius. 

Figures (6) and (7) illustrate the dependence of the 
peak pressure and the peak temperature of the bubble 
to the driving pressure and to the ambient radius, for 
the two boundary equations. Figure (6) is for the case 
in which, the ambient radius is constant (Rq = 1.5 yum), 
while the amplitude of the driving pressure is increased 
from P a — 1.4 aim to 1.8 aim. On the other hand, Fig. 
(7) shows the results for the constant amplitude of the 
driving pressure, (P a — 1-7 aim), with varying Rq from 
1.0 fim to 5.0 [im. Different values of Rq can be experi- 
mentally adapted for a specific value of P ai adjusting the 
concentration of the dissolved gas in the liquid ■ 
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FIG. 6: (a) The maximum dimensionless gas pressure versus 
different amplitudes of the pressure field according to the com- 
pressible (squares) and incompressible (filled triangles) cases, 
(b) The maximum gas temperature corresponding to the same 
data points of the case (a). The calculated points are for con- 
stant ambient radius R — 1.5 yon. Other constants and 
parameters were set as Fig.f 



FIG. 7: (a) The maximum dimensionless gas pressure versus 
different bubble ambient radii according to the compressible 
(squares) and incompressible (filled triangles) cases, (b) The 
maximum gas temperature corresponding to the same data 
points of the case (a). The calculated points are for constant 
driving pressure amplitude P a = 1.7 atm. Other constants 
and parameters were set as Fig.f 



The dependence of the maximum pressure and the 
maximum temperature of the bubble to P a is represented 
in Fig. (6). There are significant differences between 
these bubble properties resulted from the new boundary 
condition with those of the old one. The differences are 
more considerable for the higher driving pressures. More- 
over, the differences are more remarkable for the maxi- 
mum pressure than the maximum temperature. Note 
that, for P a = 1.8 atm, the maximum temperature us- 
ing the new equations is about 8 times greater than that 
of the old ones. While, the increase of the maximum 
pressure in this case is about 1400 times. These results 
indicate that, the compressibility effects should be much 
higher than what is considered in the incompressible case. 

The effects of the variation of Rq on the mentioned 
bubble properties are illustrated in Fig. (7). Similar con- 
siderable differences between the new and the old cases as 
in Fig. (6) are also present in this figure. The differences 
are more significant for the smaller ambient radii. 

Time grid resolution study of the problem shows that, 
the values of the bubble properties at the end of col- 
lapse are sensitive to the step size. But, this sensitivity 
disappears as the resolution is reduced sufficiently. Fig- 
ure (8), presents the step size dependence of the value of 
maximum gas pressure for the compressible case, where 
P a = 1.8 atm and Rq — 1.5 fim. It shows a diminishing 
step size dependence bellow 10~ 17 s. 

It should be mentioned that, although the results of 



the last figures were acquired by assumption that T = 1.4 
(the bubble content was assumed to be a diatomic gas), 
however the appearance of the new bubble behavior is 
independent of the selected value of T (and gas content). 
Figure (9) shows the evolution of the bubble temperature 
for two different states; T — 1.67 (monoatomic gas) and 
r = 1.33 (polyatomic gas). The comparison of Figs. (3) 
and (9) shows that, the new bubble behavior will be es- 
tablished even when the gas content changes. Of course, 
as Figs. (3) and (9) indicates, the configuration and the 
number of the new peaks depend on the value of T. The 
number of peaks is more for the smaller T. While, the 
values of the peak temperatures are greater for higher 
T, due to the decrease of the number of the degrees of 
freedom. 



V. DISCUSSION 

Although, the isothermal-adiabatic model used in this 
study does not account for the effects of gas dynamics, 
chemical reactions and water vapor exchange, but the 
consideration of these effects can not cover the impor- 
tance of the new term at the end of the collapse. In this 
section, the influences of these corrections on the new 
bubble behavior are argued: 

(i) Gas dynamics: Considerations of the gas dynamics ef- 
fects inside the bubble have been presented by several dif- 
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FIG. 8: Time grid resolution dependence of the maximum gas 
pressure for the compressible case, where P a = 1.8 atm and 
Ro = 1.5 iim. 



4> 4 

I- 

h- 2 





-coM.(r=i.33} 
_coM.(r=i.67) ,i 

-INCOM.(r=1.67) j| 
-INC0M.(r=1.33) ji 

/ i 
f J 



t-t^ (ps) 

FIG. 9: Dependence of the gas temperature evolution near the 
minimum radius to the value of F, according to compressible 
(thick) and incompressible (thin) cases. The dashed and solid 
curves are related to T = 1.33 and T = 1.67, respectively. 
Other constants and parameters are the same as Fig. 1. Note 
that tmin for T = 1.33 is nearly half picosecond less than that 
of T = 1.67. 



peak temperature due to the consideration of chemical 
reactions. Since most of the reactions at the collapse are 
endothermic, their influences can be considered as addi- 
tion of some extra degrees of freedom ^8|. This means 
that, the bubble actually evolves near the minimum ra- 
dius by an effective exponent, Y e ff, which is less than the 
monoatomic exponent (T — 5/3). However, according 
to our results, the appearance of the new behaviors is 
independent of the values of T. Therefore, the new bub- 
ble behaviors should be established even when chemical 
reactions are introduced in the model, 
(ii) Water vapor: Evaporation and condensation of the 
water vapor between the bubble and the liquid occur dur- 
ing the expansion and compression of the bubble. Recent 
simulations [l3l.l24l. l35l show that, a large amount of wa- 
ter vapor evaporates into the bubble during the expan- 
sion. Indeed, at maximum radius, about 90% of the bub- 
ble content is water vapor. During the collapse, the water 
vapor molecules condense to the liquid so that, near the 
minimum radius only a small fraction of a sonolumines- 
cence bubble is water vapor and the remaining is mainly 
inert gases [lj, |2J, l35| . Presence of these water vapor 
molecules inside the bubble decreases the maximum tem- 
perature due to the increase of the bubble's total degrees 
of freedom. This effect causes more decrease of the effec- 
tive exponent [l8j. Therefore, similar to our argument 
in the last paragraph, the effects of the water vapor can 
not also cover the strong effects of the compressibility 
consideration of this work. 



VI. CONCLUSIONS 



ferent approaches ill mm mill. These approaches 

show that, the overall bubble temperature and pressure 
at the collapse in gas dynamics models are at least no less 
than those of the isothermal- adiabatic model. In fact, 
in the gas dynamics models the bubble experiences more 
compression. Therefore, the pressure upon the liquid will 
be greater in the gas dynamic models. This means that, 
the effects of the new viscous term should even be more 
remarkable in the gas dynamics models, 
(ii) Chemical reactions: The temperature of the bubble 
at the end of the collapse is high enough (higher than 
10000 K) to destroy the chemical bonds of N2 and O2 
molecules of the air bubble. The chemical reactions be- 
tween the water vapor and the existing oxygen and ni- 
trogen atoms mainly produce very soluble substances in 
water (HN, NH3, and HN03), which are completely ab- 
sorbed. Therefore, a sonoluminescence air bubble mainly 
contains inert gases. This idea, which was initially pre- 
sented by Lohse et al. 33] is known as dissociation hy- 
pothesis (DH) and has been confirmed very well by ex- 
perimental reports |34j . The simulation of chemical re- 
actions at the collapse in a gas dynamics model has been 
recently presented by Storey and Szeri |24| . Their work 
shows that, a considerable decrease appears in the bubble 



The modification of the bubble dynamics equation to 
account for the viscosity of a compressible liquid was per- 
formed by deriving a new equation for the bubble bound- 
ary. This equation includes a new term, which has been 
resulted from simultaneous effects of viscosity and com- 
pressibility of the liquid. The new term is the promi- 
nent term at the end of the collapse, where the bubble 
is highly compressed. It exhibits its role by intensifying 
the strength of the collapse. The more intense the col- 
lapse is, the more significant the role of the new term 
is. Moreover, the new effects can not completely be cov- 
ered by the dissipating effects of the water vapor and the 
chemical reactions. 

The results of this work strongly indicate that, the ne- 
glect of the new term at the end of the collapse in the 
previously derived equations is not reasonable. Espe- 
cially, it is more remarkable for high amplitudes single 
bubble sonoluminescecne. It is expected that, these new 
theoretical results can be confirmed in the experiments, 
if resolution of the bubble motion measurement at the 
end of an intense collapse is less than 0.1 ns. Of course, 
a stable high amplitude sonoluminescence bubble can be 
produced, if the concentration of the dissolved gas in the 
liquid is sufficiently small |3l| . 
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